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presence of a magnetic field along the wire 
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The band structure of a quantum wire with the Rashba spin-orbit coupling develops a pseudogap 
in the presence of a magnetic field along the wire. In such a system spin mixing at the Fermi 
wavevectors —kp and kp can be different. We have investigated theoretically the collective mode 
of this system, and found that the velocity of this collective excitation depends sensitively on the 
strength of the Rashba spin-orbit interaction and magnetic field. Our result suggests that the 
strength of the spin-orbit interaction can be determined from the measurement of the velocity. 
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FIG. 1: The geometry of a quantum wire with a magnetic 
field along the wire. 



Recently active research is taking place on how to ma- 
nipulate spin properties of single electrons, and several 
semiconductor spin devices based on spin-orbit coupling 
have been proposedi^ Among these we focus on a spin 
filter^ proposed by Stfeda and Seba.* They proposed a 
spin filter combining strong Rashba spin-orbit interac- 
tion (SOI) and the magnetic field parallel to a quantum 
wire (see Fig. This system has an interesting one- 
dimensional band structure, ( see Fig. |21 ): a pseudo- 
gap is present at zero wavenumber and the orientation of 
electron spin depends on the wavevector^ For the lower 
band the electron with sufficiently negative k is mostly 
polarized in the +z direction while that of sufficiently 
positive k is mostly polarized in the -z direction. When 
the Fermi energy lies in the pseudogap substantial spin- 
mixing exists for moderate value of the Fermi energy. 
The transmission/reflection coefficients of such a wire in 
the presence of a step potential has been calculated in 



the presence of electron-electron interaction using poor 
man's renormalization group approachi^ 

The dispersion of the collective mode of quantum wires 
in the presence of a magnetic field perpendicular to the 
wire has been investigated for many years^ Recently the 
interplay of Rashba SOI and electron-electron interaction 
in quantum wires have been studied by several groupsi^ 
However, none of these studies have dealt with the case 
where the applied magnetic field is parallel to quantum 
wire in the presence of spin orbit interaction. In this pa- 
per we investigate how the collective electronic proper- 
ties may be manipulated by spin-orbit coupling. In II- VI 
semiconductors the Rashba term is expected to be larger 
than the the Dresselhaus coupling. In III-V semiconduc- 
tors, such as GaAs, the opposite is true.^ However, in 
these quantum wires the Dresselhaus term can be rather 
small under certain conditions as we argue below. The 
band structure of such quantum wires in the presence of 
a parallel magnetic field is as displayed in Fig.|2| The na- 
ture of the collective mode is unclear when the spin mix- 
ing at the Fermi wavevectors —kp and kp are different. 
We have obtained, employing bosonization method a^'^" , 
the exact dispersion relation of the collective mode of the 
lower band when the Fermi energy lies in the pseudogap. 
The dispersion relation of this mode is 

w = [vg{q)v^{q)] ^^"^q = v^q. 
V0{q) and ^^(g) are defined as follows: 

Va gV2kF 
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where uj? is the Fermi velocity and the renormalization 
factor of the strength of backscattering is 

9 = 2 ^/'^ . ^2 - (3) 
e| + (r]Rkpy 
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FIG. 2: Upper figure: Solid lines represent the lowest energy 
subband structure of the quantum wire in the absence of the 
Dresselhaus term. (Dashed lines are for zero magnetic field). 
Note that the Fermi energy lies in the pseudogap. When a 
finite value of magnetic field is present bands do anticross 
(In the figure B = 3T). The input parameters are r]R = 
2 X 10"'' eV ■ cm, m* — 0.024me. In this case the numerical 
value of g is approximately 0.7. Lower figure: The spin-up 
(solid line) and -down (dashed line) components (u^ )^ and 
(w^)^ for the lower E-{k) band. The input parameters are 
identical with the above figure. Note that (v^)'^ = 



j^gg^j^g|iiji2ji3ji4 suggest that the measurement of 
the velocity of the collective excitation vq may provide 
another way to determine the value of rjj^. This mea- 
surement can be carried out using tunneling between 
two parallel wires in the presence of an additional mag- 
netic field Bt = V X At along the y-axia^^i^^. This 
method allows one to determine the spectrum of elemen- 
tary excitations^^'^^ momenta much larger than 2kp^^. 
k [10 cm"^] This paper is organized as follows. In Sec. II we in- 

troduce our model and review the results obtained by 
Sti'eda and Seba for the non-interacting case. In Sec. Ill 
we incorporate the electron-electron interaction and ob- 
tain an effective Hamiltonian for the system. In Sec. IV 
the dispersion of collective excitation is computed based 
on the effective action obtained in Sec. III. Sec. V we 
discuss how our result for the velocity differs from the 
results of ordinary Luttingcr liquids. An experiment is 
proposed to measure 77^. 



II. MODEL FOR SINGLE PARTICLE 
HAMILTONIAN 



In our model confinement potentials are present along 
the y- and z-axis and quasi-one-dimensional motion of 
electrons is possible along the x-axis. The widths of the 
wavefunction along both the y- and z-axis are assumed 
to be negligible. The lowest subband energies along the 
y- and z- axis are denoted by Ey and E^. A magnetic 
field parallel to the quantum wire along x-axis is present 
B — —Bx. The corresponding vector potential can be 
chosen to be A = —Byz, B > 0. In our model Rashba 
electric field is applied along the y-axis (see Fig. and 
_ ^ is given by E = +Eoy, {Eq > 0). The Rashba spin-orbit 
(""fc ) • interaction^°'^^ then takes the form 



6 and (jj are the phase fields which are basically linear 
combinations of density operators Pr/l and they are de- 
fined in Eq. H33|) . Vq is the interaction matrix element, ez 
is the magnitude of Zeeman coupling, and rja is a param- 
eter characterizing Rashba SOI (see below). From the 
expression of g [Eq. Q ]we see that the velocity of this 
collective excitation depends sensitively on the Rashba 
SOI and magnetic field. This result differs from the that 
of an ordinary Luttinger liquid in that the back scatter- 
ing term V2kF is renormalized by a factor g. The physical 
origin of this factor reflects the different spin mixing of 
single particle states near the Fermi wavevectors, which 
are coupled by backscattering. 

The presence of the renormalization factor g may be 
exploited to determine the constant rjn- There is no 
simple way to calculate rjR because it depends both on 
the electric field inside the semiconductor heterostruc- 
ture and on the detailed boundary conditions at the 
interface. Instead these spin-orbit coupling constants 
were measured by electric, optical, and photoelectrical 



(4) 



where rj^ = |e|ft^i?o/4mgC^ > 0. The strength of Rashba 
SOI can be controlled by changing electric iieldm^^ 
Note that in quantum wires with electron propagating 
along the x-axis ky and fc^ must be replaced by dynam- 
ical momentum operators. The expectation value of ky, 
kz+eAz/Hc with respect to the lowest subband state wave 
function of transverse degrees of freedom {y, z) vanish by 
symmetry considerations^ 

'Hr = riiikx(Jz. (5) 

The hulk Hamiltonian of Dresselhaus SOI is given bjsi 



Wbulk.D — 7c 



^xkxiky kz) 



'^yky{k1 — k^.) + azkz{k'^ — fc^) 



(6) 



To obtain the effective Hamiltonian of quantum wire we 
have to take the average of the above bulk Hamiltonian 



3 



with respect to the ground state wave function of trans- 
verse (y, z) degrees of freedom. In our geometry the 
Rashba electric field is applied in y-direction, and the lat- 
eral confining potential enforcing quasi one-dimensional 
motion is applied in z-dircction. Clearly (kz) = since 
the subband wavcfunction along the z-axis has even par- 
ity. The subband wavcfunction along the y-axis is a real 



function and therefore the expectation value {ky 



0, 



too. But we have to note that (y) ^ since the inver- 
sion symmetry is lacking in the y-direction. The effective 
Hamiltonian for quantum wire is then 



no ^ ic<Jxt,{{kl) ~ {kl)) 



VDCTxkx 



(7) 



where r]D =lci{kl) - (kl)). ^ 

Now the one-particle Hamiltonian becomes 



we have (fc^) = m*u!o/2h. For the y-direction the con- 
stant Rashba electric field is acting so that the potential 

is linearly rising. In this case2Si (^ky) 0.8 

The condition (ky) = (fc^) is satisfied when the value of 

the electric field is given by eEoZQ = 0-^^^^^, where 

zq — y^h/'m*ujQ. For this particular value of the elec- 
tric field the Dresselhaus term is negligible compared to 
the Zeeman energy and the Rashba coupling. Note that 
E^{k) becomes an even function of k in this case. Here- 
after we assume this. If the Rashba coupling becomes 
sufficiently strong such that 



(16) 



then the energy spectrum develops a double minium at 



— Ey 



E., 



fee 

2m* 



■nnkf^z + rjokax - Ezcfx- (8) 



The Dresselhaus term can be absorbed into the Zeeman 
term Ez ~ gQiiBB/2 {go w 15 for InAs) in the following 
way. 



ez ^ Ez - rjok. 



(9) 



For the sake of completeness we include the Dresselhaus 
term in the calculation of the band structure. Later we 
will ignore it in the bosonization procedure. By the diag- 
onalization of the Hamiltonian Eq. (jHJ the energy eigen- 
values and the corresponding normalized eigenvectors are 
obtained as follows: For the lower band the eigenvalue is 
{Ey, Ez put to zero) 



2 1,2 



E-{k) = 



2m* 



rjj^k'^ 



and the eigenvector is 



where 



V(^nRk + W+4 

VRk + D 

^{mk + Dy + ei' 



Here 



D 



{mkf 



(10) 

(11) 

(12) 

(13) 
(14) 

(15) 



and represents the amplitudes for the spin to point 
in the +z and -z direction, respectively. For the upper 
band the results are given in Ref. [23 |. 

Quantum wires can be tailor made so that the quan- 
tities {ky) and (fc^) are almost equal. If we assume the 
harmonic confining potential m*ij^z^ 12 along the z-axis 



k = ±- 



1 r 

m 



* 9 

m r\ j^ 



1/2 



The energy at the minimum is given by 



E-n 



m*rij^ 



2h? 2m*ri\' 



(17) 



(18) 



In such a case -B-(O) > -Emin- In our work we assume 



that £;_(0) - £;„ 



-ez — £'min is Icss than the Fermi 



energy so that there are only two Fermi wavevectors. 



III. MODEL FOR MANY-BODY 
HAMILTONIAN 

Let and hk be the quasiparticle operators corre- 
sponding to E-{k) and E+{k) , respectively. They can 
be explicitly expressed in terms of electron operators as 
follows: 



t„,- 

fe ■ 



(19) 



When electrons are filled such that the Fermi energy is 
located in the energy gap between a and b bands at fc = 0, 
we can safely neglect the 5-type quasiparticles in the low 
energy regime. Then the Eq. (|19|l can be simplified. 



„t 



-k] "'k'^k ' '-fei "fe^fc ■ (20) 

A general electron-electron interaction in a quantum wire 
is given by 



E 



(21) 



where Vq is the interaction matrix element 



Note that 

this interaction is spin-conserving. For the long-range 
Coulomb interaction the interaction matrix element is 



2e^ 2e^ 
Vq = —Ko{\q\w) ^ — In 



1 



q\w 



for \q\w < 1. (22) 
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Kq is the modified Bessel function and w is the cutoff 
length scale which is the order of the width of the quan- 
tum wire, e is the bulk dielectric constant. For the short 
range interaction such as screened Coulomb interaction 
the matrix element Vq can be taken to be independent 
of the momentum transfer q. Projecting the Hamiltonian 
Eq. (|2HI to the a-band with the use of Eq. H2()|l . we obtain 



E 

fei ,fc2 ,q 



{ki,k2\V\ki-q, k2+q) 



where 



{ki,k2\V\ki ~q,k2 + q) 
= Vp[t{ki)C-{ki-q)][t{k2)^-ik2 



-q)] 



(23) 



(24) 



is the projected interaction matrix element in the low en- 
ergy Hilbert space. The explicit expression of eigenvector 
^_ is given by Eqs. I|11I12|I . 

At low energy, only the electron states near —kp and 
kp Fermi points need to be considered. Following the 
usual procedures of g-ology and bosonization method^ 
we can express the interaction Hamiltonian Eq. H23|) 
within g-ology scheme. Forward scattering g2 and pro- 
cess. Backscattering gi process We note further that for 
fermions of a single species (like a-quasiparticle here) gi 
process is identical with 172 process.^ In this paper a com- 
mensurate filling is not considered, so that the Umklapp 
processes (173) can be neglected. From now on we will 
call electrons with k < {k > 0) left (right) movers. 

(74 process: For instance let us assume that all four mo- 



R 



R 



[el(fci)?-(fci - q)] « [t{kF)^-{kF)] = 1. Thus the ef- 
feet of spin mixing reflected in the matrix elements 
does not play any role for 54 process. 

The contributions from the neighborhood of left Fermi 
point can be treated in the same way. The left moving 
Dirac fermion operator ip^ can be introduced similarly. 
ipLiPi) = o-ki for ki = —kp—pi. The low energy effective 
Hamiltonian describing (74 process can be read off from 
the original Hamiltonian Eq. H23|l . 



Ho 



2N 



J2 ^q\pRiq)PRi-l)+PLi-<l)PL{q)\, (25) 



k|<A 



where PR/L{q) = Y.p^'^R/L^P + q)^R/L{p) are the density 
operators of right and left moving Dirac fermions. N is 
the number of lattice sites of quantum wire. In the above 
expression the low-momentum asymptotics of Vq must be 
used. 

g2 processes: 



R 



L k. 



R 





kj-q 



kj+q 



FIG. 4; A Feynman diagrams of the (72 process. See text 
for details. There is another (72 Feynman diagram in which 
L. 



R 





kj+q 



FIG. 3; A Feynman diagram of the (74 process. All four mo- 
menta ki, k2, k2 + q, ki — q axe located near the right Fermi 
point. The dotted line indicates the matrix element Vq. See 
text for details. 



menta fci, fc2, fe-l-Q, ki — q are located near the right Fermi 
point in the following. Then one can write ki = +kp +pi 
with a condition < A < fc^- A is the momentum cut- 
off scale, within which the linearization of the a-band 
dispersion is valid. It is convenient to introduce the 
right moving Dirac fermion operator 'ip^: ipB,{Pi) = o,ki 
for ki = +kp -\- Pi. We can make following approx- 
imation if we neglect subleading contributions propor- 
tional to (k — kp) which are irrelevant at low energy: 



According to the same reason as (74 interaction one can 
make the approximations in Fig. 01 

K^(fcl)e-(fci - q)] « [t{kp)i-{kF)] = 1 

[t{k2)i-{k2+q)] ~ [t{-kp)^-{-kp)] = 1 

Again the spin-mixing effect represented by the matrix 
elements does not modify the interaction. Due to this 
one can easily read off the low energy effective Hamilto- 
nian describing 92 process from the original Hamiltonian 
Eq. H23|) using the definition of density operators PR/L{q)- 



'^32^^ E ^iPRiq)PLi-q) 



(26) 



k|<A 



gi II processes: 



In this case the matrix elements play crucial role as 
can be seen in 



[tiki)^-iki^qmUk2)^-{k2+q)] 

« (-fcF)e- i+kp)] it i+kpK- i-kp)] 



(27) 
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FIG. 5: A Feynman diagram of the giji processes. There 
exists one more diagram where _R <-> L. See text. 



Evidently the dominant momentum transfer q must be 
2kp. By changing the order of operators (thereby chang- 
ing overaU sign of interaction) and by summing over mo- 
menta one arrives at 



91 



-gV2kF / dxpR{x)pL{x), 



(28) 



where 



9 = [t{-kF)i-{+kFml{+kF)^-{-kF)] 



(29) 



4 + (rjRkFy 



As can be seen from Eq. H29|l the couphng constant g 
depend on the apphed magnetic field and the Rashba SOI 
as well as Fermi momentum. 

We observe that g2 Hamiltonian Eq. 1)26(1 and gi Hamil- 
tonian Eq. 1(281) can be combined completely. This is a 
special feature of fermions of single species. In the pres- 
ence of other degrees of freedom such as spin a backscat- 
tering term {gi.± ) appears which is not of the Luttinger 
interaction form. 



IV. BOSONIZATION AND COLLECTIVE 
EXCITATIONS 



The total effective Hamiltonian incorporating interac- 
tion is given by 



(30) 



The linearized non-interacting Hamiltonian Tio is 

= E [^FP^pi{p)i^R[p) - vfP^{[p)^l{p)\ ■ (31) 



The Hamiltonian Eq. ^ can be bosonized 
straightforwardly.^ 



n = .VFjd.[:pU.): + :plix) 
1 



2N 



X! ^9 PRii)PRi-i) + PL{-q)pL{q) 



PB{q)pL{-q)- 



(32) 



:: denotes normal ordering of operators. It is convenient 
to introduce phase fields as follows: 



{x) = ^ (I)r{x) + (t)L{x) 

[x) = \ <I^B.{x) - (t)L{x) 



(33) 



where Pr/l[x) = ^ dx4>R/L{x)- In terms of phase fields 



H^^^ I dx 
2n 



E 



V,q^ 



1 



N ^ (27r)2 



1 \^ 1 



+ (^) ^Y.^y,~9y2kM^\Q{q)o{-q)-m<i>{-q) 

+ E [2(^9 - gV2uj2)q^e{q)e{-q) 

+ gV2kFq^(biq)H~q) 
The Euclidean action is given by 



(34) 



S[9,(l>]^ I dr j dx-d^cjjdrO + H 



In matrix form the above can be written as 
1 f dujdq 



(35) 



S= — 



27r J (27r)2 

ve{q)q^ iquj 
iquj v^{q)q^ 



[9{-q,-uj) (t^{-q,-Lo)\ 



'e(q,Lo) 
<i)(q,uj) 



where 



. ^ (,^^1 V2k, \ 

ve = V9{q) = vf\1^ g- , 

\ TTVp ZTTVp ) 

.. = ..(.)=..(l + |^). (36) 

9 and are the phase fields which are basically linear 
combination of density operators PFijh ^J^d they are de- 
fined in Eq. 
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The dispersion relation of the coUective excitation can 
be obtained from the kernel of action Eq. H36|) . 



det 







After analytic continuation icu ^ lo we find 



w = [veiq)v4,{q)Y^^q = Voq 



(37) 



(38) 



Vq is the velocity of collective excitation. From Eq. H3()|) 
one can write 



TTVp 



[Vg - 9V2kj2)igV2kj2) 



1/2 



(39) 



The quantity in the bracket of Eq. H39() represents the 
renormalization effect due to electron-electron interac- 
tion. The velocity of collective excitation can be con- 
trolled by band filling, Rashba SOI, and magnetic field 
through dependence on vp and g. Let us estimate 
the magnitude of the correction terms. In vg^q) the 

backscattering term, g , is a factor g/2 smaller than 
the forward term, In v^{q) the correction term 

gV2kF /'2ttvp ^ O.I5 for the width of the quantum wire 
w ^ lOOA and 2fci? w 1 x lO^cm"^. We also note that for 
the screened short range Coulomb interaction the inter- 
action matrix element Vq is almost independent of mo- 
mentum transfer q, and the backscattering term plays an 
equally important role as forward the scattering. 



V. DISCUSSIONS AND SUMMARY 

It is instructive to compare this result with the veloc- 
ities of phase fields of collective excitation of ordinary 
Luttingcr liquids. For spinless fermions it is given by 

ve{q) = 1 + - 7y , 

\ TTVp ZlTVp ) 



v^{q) :^vp[l + 



2nvp I 



(40) 



In Eq. (|2Jl this corresponds to g = 1, which implies 
absence of spin-orbit coupling and one type of spin, either 



up or down. For the Luttinger liquids of spinful fermions 
the velocity of charge mode is given by 



, 2V„ 
ve„ {q) = vp\l + 



2kF 



TTVp 2t:vp J ' 
...(.) = + (41) 



6„ and 



p are the phase fields in the charge sector. 
The spinful velocity is recovered with the replacement 
Vq 2Vq and g = 1 in Eq. (0). The velocity of the spin 
mode is 



^e, (q) =vp [1 



v. 



2kf 



(<?) = I'-F 1 + 



27rvp J 



(42) 



2nvp J 

Os and (ps are the phase fields in the spin sector. This 
corresponds to Vq = Q and g = 1 in Eq. 

The dispersion relation of the collective mode may be 
measured by adding another quantum wire parallel to the 
original wire in the presence of a second magnetic field Bt 
along the y-axis. When the first wire is located at z = 
and the second wire at z = the single particle energy 



h^(k~koY 
2m 



where 



dispersion of the second wire is E{k) 
fco — eBtZQ/hc, m is the electron mass in the second wire, 
and the Landau gauge A( = {zBt, 0,0) is used. Wave- 
number selectivity due to momentum-resolved tunneling 
between them, E-{k) = E{k), allows a mapping of the 
dispersionpiSiiSiiS. Even in the presence of electron inter- 
actions this technique allows direct measurement of the 
collective excitation spectrumiiLi^ 
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